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^ Abstract 

o 

CN Simulation of the CERN crystal-extraction exper- 
^iments at the SPS [1], including the real crystal 
Q geometry and SPS parameters, as well as multi- 
pie turns and passes, has been performed. Two 
CN| crystals with different geometry have been tried. 

yVe argue that in both cases the extraction ex- 
;> periment suffered from the crystal edge imperfec- 
r"- tion. For the experimental conditions reproduced, 
^the efficiency of the first ("twisted") crystal was 
O found to be 12-18% (peak) in the angular range of 
140-260 /xrad (FWHM), depending on the vertical 
O beam position at the crystal location. For the sec- 
^ ond ("U-shaped") crystal the peak efficiency was 
^much the same, ~19 %, while the angular scan 
QhFWHM reduced to 70 /xrad. These results, as 
^ c-| well as the beam profiles, are in agreement with 
^ measurements. Finally we suggest some ways to 
advance this experiment. 

n\ 

1 Introduction 

The series of crystal extraction experiments is be- 
ing performed at CERNjl]], studying the feasi- 
bility of extracting protons from the halo of the 
beam circulating in SPS by a bent silicon crystal. 
These studies aim to model the beam parameters 
expected at the future Large Hadron Collider, and 
have in view possible application of the technique 
for a beam extraction from a multi-TeV machine. 
Such an extrapolation requires full understanding 
of the crystal extraction process. The present pa- 
per investigates these experiments by means of a 
computer simulation. Our aim is not only to com- 
pare theory vs measurement, but also to discover 
the factors apparently disturbing the crystal ex- 
traction process in M. 



The physics of channeling in a bent crystal is 
well established [0, |3|. Beam bending by a crys- 
tal is due to the trapping of some particles in 
the potential well formed by the field of atomic 
planes, where the particles then follow the di- 
rection of (are channeling in) the atomic planes. 
The channeling effect survives in a bent crystal 
until the ratio of the beam momentum p to the 
bending radius R becomes as high as the maxi- 
mal field gradient (~6 GeV/cm in silicon). The 
crystal bend reduces the phase space available for 
channeling, thus decreasing the fraction of par- 
ticles channeled. The scattering processes may 
cause the trapped particle to come to a free state 
(dechanneling). 

In the present simulation we tracked protons 
through the curved crystal lattice with small, ~1 
(Urn, steps applying the Monte Carlo code CATCH 
[|J. This code uses the Lindhard's continuous- 
potential approach to the field of atomic planes, 
and takes the processes of both single and mul- 
tiple scattering on electrons and nuclei into ac- 
count. We assumed the crystal to have a perfect 
lattice. 

The simulation of extraction was performed 
with parameters matching those of the SPS ex- 
periment: an= 2.07, (3h= 90 m, Qh— 0.635, 
a v = -0.734, (3 V = 24.4 m, Q v = 0.583. The 
vertical normalized emittance ey was normally 
10 7rmm-mrad. Protons were tracked both in the 
crystal and in the accelerator for many subsequent 
passes and turns before eventual loss either at 
the aperture or in interaction with crystal nuclei. 
Tracking in the SPS used linear, 4 by 4 transfer 
matrices. The crystal was located horizontally at 
10 mm from the beam axis. 

The beam was excited by horizontal white noise 
at the opposite side of the machine, with kicks of 



0.003 /irad r.m.s. value per turn. The expected 
impact parameters b and angles b' were studied 
earlier ]5|. The b' value there was found to be 
<C ip p . The b value should be ~0.7/im, i.e. quite 
close to the surface. 

The crystal surface quality setting was rather 
conservative: miscut angle 200 /irad, surface non- 
flatness 1 /im, plus 1 yum thick amorphous layer 
superposed. This defines some 'septum width' of 
a few /im value. The protons tracking near the 
surface have taken all these details (as well as the 
bent surface) into account. We consider here two 
options. The first, with impact parameter below 

1 /im and surface setting as described above, ex- 
cludes the possibility of channelling in the first 
pass through the crystal. This is compared to the 
second option, when the crystal surface is assumed 
perfect, i.e. with zero septum width. 

2 "Twisted" crystal 

The first extraction experiments [ 1 ] were per- 
formed with a " Serpukhov- type" bending device, 
and have given results more sophisticated than ex- 
pected: the extraction efficiency was rather small, 
around 10%; extraction occurred in the angular 
range of 200 /irad FWHM , much wider than the 
beam divergence or critical channelling angle (14 
/irad); two peaks appeared both in the horizon- 
tal and vertical profile of the extracted beam at 
the crystal orientation far away from the best one; 
the angular range of extraction depended on the 
vertical beam position. 

In the laser measurements it was found that the 
bent crystal is deformed vertically, called 'twist' 
in the following. As a result, the planar direction 
at the entrance to the crystal was a function of 
the vertical position. The measured crystal prop- 
erties were used as input to the simulations. The 
total bending found for this (3 cm long, 1.5 mm 
thick) crystal was 8.5mrad. In the experiment, 
different vertical beam positions (bump) and dif- 
ferent vertical emittances ey have been tried in 
order to investigate the twist effect. Our simula- 
tion followed the program of the SPS experiment. 
Table 1 summarizes the extraction efficiencies, the 
widths of the angular scan and of the extracted- 
beam profiles (observed 20 m downstream) ob- 



Table 1: Summary of simulation for ey=107r and 

1.57T. 
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tained in the simulation for ey = 107r and 1.57T 
mm-mrad. The extraction efficiency function 
of tilt angle is shown in Fig. 1 both for nominal 
position (•) and bump of 2 mm (o). The increase 
of ey to 20 it mm-mrad has slightly reduced the 
efficiencies. This is summarized in Table 2, where 
the peak efficiency is plotted as a function of ey 
and the vertical position of the beam at crystal. 

In the SPS experiment, the nominal vertical po- 
sition was ~lmm away from the 'twist middle'. 
The ey typical value was ~107r mm-mrad. There- 
fore the number of 16% for the peak efficiency 
(and 170 /irad for the angular scan FWHM) seems 
the most relevant one to compare with the exper- 
imental results, namely 10±1.7% and 200 /irad. 
The systematic error caused by some uncertain- 
ties in the vertical position, ey value, and twist 
detail, can be estimated from Table 2 few 
per cent for the peak efficiency. We may conclude 
that the perfect-surface simulation has often pre- 



Table 2: The peak efficiency (%) plotted vs ey 
and the vertical beam position (bump). 
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Figure 1: Angular scans of efficiency for twisted 
crystal. Nominal position (•) and bump of 2 mm 

(o) 



dieted narrow high peaks, which have not been 
observed. The imperfect-surface option roughly 
followed the experimental observation. Both the 
angular scan FWHM and the efficiency value de- 
pend on the bump in the same way as it was ob- 
served at the SPS. The extracted beam vertical 
width tends to decrease with vertical displace- 
ment, while the horizontal width has no trend but 
sizeable fluctuations. Both features seem to agree 
with experimental observation. There is quanti- 
tative agreement with the measurement, but the 
price for this agreement is the assumption of crys- 
tal edge imperfection. The twist effect was minor 
for the peak efficiency, and cannot alone explain 
the large angular width of the extraction. 

3 "U-shaped" crystal 

After the above analysis, a prediction was made in 
for the new, "U-shaped" crystal H] which has 
no twist. This crystal is 4 cm long with 8.2 mrad 
bending. Fig. 2 shows the angular scan simulated 
for this crystal with edge imperfection, which is 
in good agreement with the first experimental 
results [0]. As compared to the "twisted" crys- 
tal, there was expected no sizeable change of effi- 
ciency, in accord with the first observation. Such 
an agreement indicates that the new crystal also 
has some edge imperfection. For an ideal crystal 
and a very parallel incident beam the simulation 
gives the peak efficiency ~50% and the narrow 



angular scan (25/zrad FWHM). 

4 Experiment Advancement 

An essential advancement of the crystal extrac- 
tion experiment would be the measurement of the 
first-pass contribution to efficiency. This is quite 
important for understanding both of the crystal- 
edge work (measuring "septum width" ) and of the 
crystal co-work with the other elements of accel- 
erator (multiple passes/turns). One can distin- 
guish the first pass from the secondary ones ow- 
ing to their difference in time (first comes first) 
and space (secondary pass comes deeper in the 
crystal). The time-synchronyzing ("kick mode") 
has difficulty caused by unavoidable mixing of the 
"first-pass" protons from later turns. For a spa- 
tial separation one way is a thin amorphous layer 
superposed on the crystal surface to suppress the 
first pass (W.Herr). Here are some other ways: 

Microscope. Applying a skew cut on the exit 
face of a crystal, one can make the bending angle 
to depend on the impact parameter b (a princi- 
ple of any microscope). The bent beam profile 
downstream the crystal would be then an ampli- 
fied profile of this beam at the crystal face. 

Energy loss The 5E/5z deposited in the 
crystal may indicate whether the proton is ex- 
tracted with the first pass or with the secondary 
one, due to the escape of (^-electrons near the sur- 
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Figure 2: Prediction (®, left scale) and measure- 
ment (*, right scale) for the U-shaped crystal. 




Figure 3: Extraction efficiency vs crystal length. 
For edge imperfection (•) and ideal crystal (o). 



face (depending on b). And vice versa, selecting 
'first-pass' in some other way one may observe an 
interesting SE/Sz spectrum (fwhm^O). 

Masking [6[]. One can scrape the circulating 
beam by a collimator during the crystal extrac- 
tion, thus suppressing the secondary passes. 

Tune resonance [Q. It was shown that for 
an imperfect edge of a crystal the extraction effi- 
ciency drops near the tune resonance (for the same 
given parameters of incident particles). The reso- 
nance both dip and width depend on (are roughly 
proportional to) the septum width. 

Crystal rotation in "FNAL geometry" 
[0. In the different geometry of extraction Jl0|] 
there is an extra degree of freedom. The crystal 
rotation normal to bending plane results in a con- 
trolled change (with steps of ~0.02 //m) of septum 
width. 

Fig. 3 shows the extraction efficiency as func- 
tion of L, simulated for an untwisted crystal with 
constant curvature. A short (~1 cm here) length 
could be a " cure" to any imperfection of a crystal, 
as the protons make more passes. 
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